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Abstract 

A set of moderately deformed s — d shell nuclei is employed for 
testing the reliability of the nuclear ground state wave functions which 
are obtained in the context of a BCS approach and offer a simultaneous 
consideration of deformation and pairing correlations effects. In this 
method, the mean field is assumed to be an axially symmetric Woods- 
Saxon potential and the effective two-body interaction is a monopole 
pairing force. As quantities of main interest we have chosen the nuclear 
form factors, the occupancies of the active (surface) orbits and the 
Fermi sea depletion, which provide quite good tests for microscopic 
descriptions of nuclei within many body theories. For our comparisons 
with results emerging from other similar methods, an axially deformed 
harmonic oscillator field is also utilized. 

PACS : 21.60.Jz, 23.40.Hc, 27.30+t, 21.10 Ft 
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1 Introduction. 



In recent nuclear structure calculations, the reliability of various predictions 
associated with phenomena of several prominent electroweak processes {j3- 
decay modes, lepton-capture by nuclei, neutrino-nucleus reactions, etc.) re- 
lies upon the accurate description of the nuclear many-body wave functions 
Among the nuclear structure aspects which can quite strongly influ- 
ence such predictions, the consideration of the deformation and the 
nucleon-nucleon correlations effects [3 |HJ 13 HH] (wherever they might be ap- 
preciable) is of significant importance. As it is well known, the ground state 
nucleon-nucleon correlations [7] are of particular significance to be reliably 
incorporated especially when investigating processes for which the branching 
ratios or decay (event) rates are highly suppressed and, therefore, the corre- 
sponding transition matrix elements are very sensitive [3EI- Moreover, the 
nuclear deformation, e.g. for axially symmetric nuclear systems, should nec- 
essarily be considered by applying appropriate treatments and using reliable 
methods jSlilElini- 

The main purpose of the present work is to address the effects of deforma- 
tion and pairing correlations on a few nuclear observables of the ground state 
for moderately deformed s-d shell nuclei in the spirit of approaches based on 
BCS quasi-particles. The properties of main interest are the ground state 
charge form factor Fch{<i) [HI 1121 USl UH UH); as well as the occupancy (va- 
cancy) numbers of several orbitals around the Fermi level fOl 1^ ^1 ^1 120] • 
Such quantities (also those referred to low-lying excitations like transition 
densities, transition form factors, etc.) constitute very common tests in 
nuclear structure theory for the many-body model wave functions ITB] . 
Experimentally, the above properties are known with very high accuracy 
(~ 1%), e.g. from experiments of elastic scattering of electrons, protons, 
etc. off target nuclei [HI [121 CSl El ■ On the theoretical side, average poten- 
tials (spherical Woods-Saxon or harmonic oscillator, etc.) employed in BCS 
theories and nucleon orbitals, both of which are directly related to gross nu- 
clear features of the ground states, have been quite effective in describing 
the dynamics of many mainly spherical nuclei |TmiT71l^ l21j. However, even 
though the direct calculation of the charge form factor (and other gross nu- 
clear properties) through a nuclear model ground state (g.s.) wave function 
is a very extensively discussed issue, the reproducibility of the scattering data 
in the high momentum transfer region {q > 2 — 3 fm~^) with the experimental 
accuracy has not yet been achieved (in the context of any nuclear model) for 
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most of the nuclear systems. 

The existence of nuclear deformation, in general, is closely related to the 
non-spherical components of the proton (neutron) density distribution and 
the corresponding proton (neutron) nuclear form factor whereas the nucleon- 
nucleon correlations are, to a large extent, related to the nuclear surface 
texture, the main characteristics of which are the occupation probabilities 
of the active orbits. These quantities are studied in single-nucleon transfer 
reactions [TS], but often the results differ substantially from the predictions 
of independent particle approximations and other methods, due to the fact 
that, in many of these models, the correlations and deformation phenomena 
are not (or may not correctly [22]) taken into account [7]. Hence, details of 
the nuclear surface and also bulk phenomena for many nuclei may have not 
correctly been understood in terms of many existing models. 

The present method employs a BCS ground state built out of a deformed 
WS basis wave functions. The required two-body interaction is taken to have 
a schematic-type, the monopole part of which is scaled phenomenologically 
(separately for protons and neutrons) so as to be consistent with the observed 
proton and neutron separation energies within the BCS procedure I2S]- 
At first, we study the nuclear form factors for the set of deformed nuclei 
^2Mg, i2Mg, llSi and fgS (s-d shell nuclei), by considering them as systems 
having axially symmetric f3 (quadrupole) deformation only f]M |^. They 
are rather well deformed isotopes with the following values of the parameter 
/3: PC^Mg) = 0.416 1^ (0.374 |^), p^^Mg) = 0.296 [25^ (-0.310 j211), 
/?(28si) = -0.328 [2S] (-0.478 ^) and (3{^^S) = 0.186 ^ (0.000 (211). In 
the present work we apply also values of the deformation parameter f5 which 
originated from experimental values of the quadrupole moment, taken from 
Ref. p6J 

We next proceed with the study of the effects of deformation and pairing 
interactions on the surface features of these nuclei, namely, on the smearing 
of their Fermi surface which is very sensitive to the orbit occupancies of the 
active space. In our method, deformation effects and pairing correlations for 
proton-proton and neutron- neutron interactions may simultaneously and in a 
quite reliable manner be considered which is in accord with the experimental 
indications of Ref. jT3] showing that these effects must be treated as much 
as possible on equal footing. Finally, for the sake of comparison with the 
results of other similar methods, we study the same properties for the above 
isotopes, but now by expanding the required deformed Woods-Saxon wave 
functions into a deformed harmonic oscillator (DHO) basis. 
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We remark that, spherical WS potentials parameterized for the well depth, 
skin thickness, radius and spin-orbit constants, etc. as in Ref. j27j, have been 
used by us in the past for microscopic random phase approximation (RPA) 
calculations of decay (event) rates in electroweak processes (/5/3-decay 
flavor violating processes [SI 1211 , ^^^O ^^r various nuclei. Several of these 
calculations must be extended to non-spherical systems or must be corrected 
to include deformation and/or correlations effects. The present method offers 
some advantages for such purposes |2Sl 123 • 

In the remainder of the paper, we first describe the main ingredients of 
our method (Sect. 2), then we present the results obtained and compare 
them with existing experimental data (Sect. 3), and finally (Sect. 4) we 
summarize the main conclusions extracted. 



2 Brief description of the formalism 

2.1 Axially symmetric Harmonic Oscillator and Wodds- 
Saxon bases 

For the description of axially symmetric nuclear systems, where the cylin- 
drical coordinates (r, 0) are appropriate, a deformed harmonic oscillator 
(DHO) potential has been primarily used jSU]. The eigenfunctions of this 
unharmonic oscillator are described in terms of two frequencies: a;^, in the 
direction of the z-axis, and cj^, for the motion perpendicular to the symmetry 
z-axis. Then, the space-spin DHO single particle wave functions are written 
as [SOI 

^np,n,,A,s(r,Cr) = *„p,„^,a,e(P,2:,0,O-) = V'i'^'(p) i^nA^) V^a(0) X{^) (1) 

where S and A are the projections on the z-axis of the spin and orbital 
angular momentum, respectively, x(^) stands for the spin wave function 
and the components of the spatial parts are defined as follows: 

(i) The radial part is usually written in terms of a dimensionless coordi- 
nate r] defined by 

.) = /^i^.^ flx = (^)', (2) 

(m is the nucleon mass) as 

<'(p)=<'(^)^^^e-i<l(,), (3) 
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where 

\np+\K\)\' ' 

is a normalization factor and LI^' (r^) are the associated Laguerre polynomials. 

(ii) The z-dependent part is similarly written in terms of a dimensionless 
variable ^ 

i = R.z, R. = {^)K (5) 



as 



^nM=NnA'^t^e-'TH^^{0 (6) 



where again 



= (v^2-nj)-^ (7) 

is a normalization factor and HnXC) ^-^e the Hermite polynomials, 
(iii) the 0-dependent part is given by 

M<P) = ^e^^^ (8) 
v/vr 



which is normalized to unity. 

The coordinate scaling parameters and Rz-, which have dimensions of 
inverse length, characterize the motion in the perpendicular and z-direction, 
respectively. Obviously, the level of anisotropy of the field is represented by 
the difference between the two frequencies, uj_^_ and cj^, which are written in 
terms of the usual quadrupole deformation parameter (3 as [20] 



^o(/3) (1-^/3), 0;^ = ^o(/5) (1 + 



where ujq is weakly dependent on j3 (enough to conserve the nuclear volume). 
From the latter two equations we obtain 

The parameter is defined so that > correspond to the so-called prolate 
shapes and /? < correspond to oblate shapes. 

The eigenfunctions of the Hamiltonian involving an axially symmetric 
Woods-Saxon potential have a good quantum number K, the projection of 
the total angular momentum, J = L + S, on the symmetry z-axis, i.e. K = 
A + S (due to the assumed symmetry they are eigenfunctions of the Jz 
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operator). These eigenfunctions can be expanded into basis functions like 
those of Eq. IQ, namely, into the eigenfunctions of a deformed harmonic 
oscillator potential. In the intermediate region of the parameter (3 we are 
interested in, 0.15 < /3 < 0.50, one could employ either Eq. (|T)) or total 
angular momentum coupled bases like the 'if{N,nz,A,K) (see Ref. jSU] P- 
118) which we apply in this work (here the principal quantum number N is 
defined as = + 2np+ | A |). In this case the expansion of the WS basis 
takes the form 

MK) = Yl «.(iV,n„A)^(A^,n„A,S = l/2,ir) 

+ «:(Ar',<,A')vI/(Ar',<,A' = A + l,S' = -l/2,K) (9) 

where a, and a'^ are expansion coefficients obtained from the diagonalization 
procedure. The time-reversed partners of the states have the form 

^,{K) = Y a,{N,n,,A)^{N,n,,-A,J: = -l/2,-K) 

N,n^,A 

- Y a:(iV',<,A')*(^^',<,A' = -A-l,S' = l/2,-ir|10) 

Using Eqs. Q and (fTUI). we can write the expression for the point proton 
density distribution in an axially symmetric nucleus as 

P(r) = YK[^:iK)MK) + ^:{K)MK)] (11) 

i 

where Vp^ stands for the occupation probability of the z*^ proton state (in 
the deformed basis), determined from the solutions of the BCS equations. 
The proton density distribution, p(r), Eq. (fTTj). is related to the proton form 
factor -F(q) of the nucleus in question by the well known expression (Born 
approximation) 

F(q) = J p(r)e^'i'"dr (12) 

By choosing q in the direction of z we have, iqr = iqr cos 6 = iqz. Under this 
assumption, the nuclear form factor, Eq. (fT^ . can be written in the form 

Fil) = 2E I E E c^^iN,nz,A) «,(A^',<,A) /(n„<,g) ] 

i \N,n^,AN',n'^ / 

(13) 
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where 

/(n„ K,q)= I C {z)e'''dz. (14) 



The first part of our calculations in the present work (see Sect. 3) relies upon 
Eqs. (USD and dUD. 

2.2 The BCS model 

The ground state of even-even deformed nuclei is determined by the deformed 
pairing Hamiltonian, which includes monopole proton-proton and neutron- 
neutron pairing interactions as jSJ EI] 

S (7 S (T 

s,s' s,s' 

where e^^ and are the un-renormalized single particle energies for protons 
and neutrons, respectively. \p (A„) is the proton (neutron) Fermi energy and 
the operator Sl^^, creates pairs in time reversed orbits written as 

spp psa psa ' ^ snn Z-^ nsa nsa • V 

cr cr 

Here, c\^^ and c^^a stand for the creation and annihilation operators of a 
particle (t = p and t = n denote proton and neutron, respectively) in 
the axially symmetric Woods-Saxon potential, a is the sign of the angular 
momentum projection K, namely, a = ±1. We note that the intrinsic states 
are twofold degenerate, since the states with K and —K have the same energy 
as a consequence of the time reversal invariance (the symbol ~ indicates 
time reversed states). The pairing interaction strengths Gpp and which 
characterize the associated monopole interaction for proton or neutron pairs, 
respectively, are determined by fitting the theoretical pairing gaps to the 
average empirical ones parameterized as, App = Ann = 12.84/ A^^'^ MeV fl^ . 
where A represent the mass number of the isotope in question. 

3 Results and Discussion 

The first part of the present calculations refers to the proton form factors of 
the deformed nuclei i|Mg, f^Mg, HSi and ^gS carried out as follows. At first. 
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in the spirit of the independent particle approximation (Slater determinant 
calculations), we used a deformed Woods-Saxon field parameterized as in Ref. 
P7j (in this step there is no correlations effect). In the second step, the pairing 
correlations are inserted in our calculation by means of the BCS theory, 
namely through the monopole pairing interaction, a schematic type nucleon- 
nucleon force. The strength of this pairing force, is determined separately for 
protons and neutrons as in Ref. [S] (see Sect. 2.2). In the third step of the 
calculational procedure, for comparison, we evaluated the same proton form 
factors utilizing now a deformed harmonic oscillator basis in the BCS ground 
state, as discussed in Sect. 2. Thus, the latter results include correlations 
of the same kind as those of the second step. For all the cases studied, 
we adjusted the parameter of Eq. © so as, for each isotope, the first 
diffraction minimum of the nuclear form factor to appear at the experimental 
position (value of the momentum transfer q). For the deformation parameter 
(3, throughout this part of work we used the more recent calculations of 
Ref. which are, generally, in agreement with the predictions of the 

macroscopic-microscopic model of Ref. [21] . 

The results obtained as described above, are presented in Fig. 1, where 
we use the notation WS to denote the deformed (axially symmetric) Woods- 
Saxon results, the notation WS-I-BCS to represent the results containing in 
addition pairing correlations, and the notation DHO+BCS for the results 
which consider the deformation and pairing-correlations effects through the 
harmonic oscillator model. The main conclusions stemming out of the above 
calculations are the following. As can be seen from Fig. 1, the incorporation 
of pairing correlations (second step) via the BCS procedure, affects apprecia- 
bly the shape of the charge form factors (compare the plots corresponding to 
WS and WS+BCS results), especially for q> 2 fm~^, which is a well known 
outcome. Furthermore, the simultaneous consideration of deformation and 
pairing correlations in the BCS ground state with Woods-Saxon potential, 
reproduces the experimental form factor data (and charge density distribu- 
tion) much better than the deformed harmonic oscillator model (compare 
results labeled WS+BCS and DHO+BCS). As it is evident from Fig. 1, 
nearly all the diffraction minima of the results labeled WS+BCS are pro- 
duced quite close to their experimental momentum transfer locations (also 
some additional minima not observed by experiments up to now are predicted 
for values of g > 3 — 4 fm~^) [201 • -^^^ high momentum transfer (close to the 
nuclear center of the charge density distribution) the reproducibility becomes 
worse due to the presence of other sorts of correlations not included in our 
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results (short range nucleon-nucleon correlations, etc.) and other kinds of 
effects (meson-exchange corrections, etc.). 

We investigate also the effect of the deformation on charge form factors 
as well the occupancy numbers by using various values of the deformation 
parameter (3 which originated from Ref. [20] (see also Table 1). More specif- 
ically the parameter j3 calculated through the relation 



where Qq is the intrinsic quadrupole moment (in barn), is the charge radius 
(in fm) and Z the number of the proton. Qq calculated from the experimental 
quadrupole moment Q via the relation 



The relative values are exhibited in Table 1. 

Fig. 2 focuses on the nucleus ^"^Mg and ^^S showing the calculated nuclear 
form factors for various values of the deformation parameter (3 where taken 
for Table 1. We see that, the shape of the form factor is sensitive to changes 
of the quadrupole deformation. It is worth noting that, in general, the values 
of the parameter (3 affect twofold the form factor results. First, they define 
the set of the single particle WS wave functions for Eq. (fT^ and energies 

of Eq. (fT3j) and, second, through the single-particle energies e^^ and the 
subsequently given pairing strengths, they determine the values of the orbit 
occupancies entering Eq. (fT^. 

In the next part of our work, for each isotope studied we investigated the 
influence of pairing correlations and deformation on the details of the nuclear 
surface indicated by the occupancies of the active orbits and the Fermi sea 
depletion (the sum of the valence vacancy numbers, U^^). In Fig. 3, the 
proton-orbit occupancies versus the single particle energies are illustrated 
(cases denoted as WS and WS-I-BCS) and in Table 2 (last column) we report 
the Fermi sea depletions estimated for the case WS-I-BCS of Fig. 3. From 
these results we conclude that, the largest depletion corresponds to ^^Mg 
nucleus (the most deformed isotope of our set) and the smallest one to ^^S 
(the least deformed isotope of our set) which might be accidental. Concerning 
the two Mg isotopes, evidently, the depletion of the less deformed ^^Mg is 
appreciably smaller than that of ^^Mg showing the quite strong effect of 




(17) 
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Table 1: The experimental values of the quadrupole moment Q and the 
corresponding values of the deformation parameter (3 taken from Ref. |26j . 



Nucleus Z Tc (fm) Q (barn) [3 

24Mg 12 3.126 -0.16 ^ -0.29 0.38 ^ 0.69 

26Mg 12 3.083 -0.11 ^ -0.21 0.27 ^ 0.51 

28Si 14 3.177 +0.16 +0.18 -0.31 ^ -0.35 

^^S 16 3.282 -0.12 -0.18 0.19 0.29 



deformation on the smearing of the nuclear surface. This result agrees with 
the "hardening" of the valence-proton distribution found in Ref. [ISl which 
is due to the addition of two extra neutrons. 

In Table 2, for comparison, we quote the shell model Fermi sea deple- 
tions of Refs. jini E] discussed in conjunction with elastic scattering data 
that provide information about differences in occupations of individual orbits 
between the pairs of neighboring even-even nuclei (^"^Mg, ^^Mg) and (^^Mg, 
^^Si). In this Table we also include the shell model Fermi sea depletion cal- 
culations of Ref. [HHl and the values extracted from the experimental data 
of Refs. inn ES] • The authors of Ref. [^3] , using two well-tested effective 
two-body interactions [311 EZ| and treating the above nuclei as spherically 
symmetric systems, found drastically smaller depletions compared to ours 
and appreciably smaller than the experimental ones and those of Ref. P^ . 
From Table 2 one can infer that, the experimental Fermi sea depletions of 
Ref. [34j 35j are overestimated by about 40-130 % (depending on the nuclear 
system) from our method, and they are underestimated by about 10-50 % 
from the latter shell model calculations. The calculations of Ref. [13], in 
general, agree quite well with the experimental results. 

The disagreement between these experimental data and our theory might 
be related to the factors discussed below, after making the following remarks: 
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Table 2: Fermi sea depletions for the deformed nuclear isotopes ^^Mg, ^^Mg, 
^^Si and ^^S evaluated by the present work (last column) and previous spher- 
ical shell model calculations. The experimental results are from the following 
Refs. SM I IMl, SM II HH, SM III mi, Exp. I |M1 and Exp. II This 
result comes from spherically symmetric Hartree-Fock calculations discussed 
in Ref. which include the O/7/2 orbital in the model space. 



Nucleus 


(MeV) 


SM I 


SM II 


SM III 


Exp. I 


Exp. II 


This Work 


24Mg 


2.62 


3.6 


4.0 


9.2 


6.7 


7.3 


16.8 


26Mg 


2.52 


4.3 


3.9 


6.7 


5.8 


6.0 


12.6 


28gi 


2.43 


6.9 


6.8 


9.8 


10.7 


10.7 


15.3 




2.27 


5.2 


4.8 


5.8* 


5.0 


5.6 


11.1 



In our BCS treatment all orbits involved in our rich model space are consid- 
ered as active, whereas the shell model results of Ref. ^1 ES] originate from 
the very restricted model space 0d5/2lsi/20d^/2 (both interactions used in [32] 
give rather similar orbit occupancies even though they differ in their predic- 
tions of occupancy dependent single particle energies, which means that for 
each isotope, the corresponding to our Fig. 3 pictures given by the latter two 
shell model calculations might be substantially different compared to each 
other). The parameterization of the deformed WS potential assumed in this 
paper, had originally been proposed for spherical nuclei (ranging from g^O 
to g2^Pb). The confidence level of this parameterization (and subsequently 
of our results for deformed nuclei) relies upon its isospin-dependence charac- 
teristic which has indirectly been checked from the reliability of descriptions 
done (under the same assumptions with the present work) in Ref. [3 EHl 120] • 
In order to elucidate further the relationship between the deformation 
parameter /3 and the depletion of the valence shells, we have carried out 
a set of calculations with different values of f3 (originated from Ref. ^Ei, 
see also Table 1) keeping fixed the pairing gap for the proton-proton and 
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neutron-neutron channel, App = A„„, in the BCS procedure (restricting in 
this way as much as possible the variation of pairing correlations effects). 
We have chosen the nuclei ^"^Mg and ^^S and the results are displayed in Fig. 
4. It is evident from this figure, that the Fermi sea depletion decreseased 
slightly as the deformation effects tend to be canceled out and vice versa. 
More precisely the depletion ranging in the case of ^^Mg from 17% (/3 = 0.38) 
to 17.5% {(3 = 0.69) and in case of ^^S ranging from 11.11% {/3 = 0.19) to 
11.31% {(3 = 0.29). This is an interesting result and may be accounted for as 
follows. By keeping the pairing interaction rather constant, a decrease of the 
deformation increases the orbit vacancies t/pj and subsequently causes more 
effective participation of the nuclear core to the surface texture. Moreover, 
since for the latter results (and those of Table 2) we have not considered other 
sorts of correlations except pairing ones, the smearing of the Fermi surface 
is artificially interpreted by increasing the vacancies of the hole states (nu- 
clear core), namely the Fermi surface depletion. In the region of deformation 
parameter, 0.2 < (3 < 0.7, where our method is applicable, the Fermi sea de- 
pletion seems, for this reason, to be overestimated. At this point we should 
mention that, in the recent literature JHI authors assert that the depopula- 
tion of states below the Fermi energy is much more pronounced than it was 
believed in the past, since the quantitative information we have today from 
both theory and experiment support occupation probabilities in finite nuclei 
of the order of 75 % only. Our present results underly this result and it is 
expected from future measurements on specific nuclei with higher precision 
and covering larger range of momentum transfer to shed more light on the 
situation. 

Before closing we note that, in our calculations, effects like finite pro- 
ton size, center-of-mass corrections, etc., have not been taken into account. 
Instead, we restricted ourselves to the evaluation of the point nucleons (pro- 
tons) distribution, which is equivalent to ignoring their internal substructure. 
This approximation has extensively been used in nuclear structure methods 
based on BCS ground state, because, the shape of the charge form factor is 
mainly determined by the shape of the point-proton density. In the present 
work we have also stayed within the assumptions of the above approxima- 
tion encouraged for it by the rather well description of decay processes like 
the double beta decay |29j performed recently by using wave functions con- 
structed as described in Sect. 2. As it is well known, the nuclear structure 
details affect considerably the corresponding double beta decay rates since 
its transition matrix elements are fairly sensitive. 
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4 Summary and Conclusions 



In the present paper, at first, we formulated on the basis of BCS quasi- 
particles a simultaneous consideration of deformation and pairing correla- 
tions phenomena. The mean field description is provided in terms of an 
axially symmetric Woods-Saxon potential whereas for the two-body interac- 
tion we assumed schematic forces. Next, wc tested this method by studying 
gross ground state properties like the nuclear form factor, the occupancy 
(vacancy) numbers of the active orbitals, and the corresponding depletion of 
the Fermi sea for the set of deformed nuclei, 12 Mg, 12 Mg, ft Si and fgS. As it 
is well known from nuclear structure calculations, the focus on such primary 
gross properties, which is by itself a continuously important subject of study, 
provides very good tests for microscopic descriptions of the nuclear dynamics 
within many-body theories. This is encouraged by the great resource of data 
for these properties provided mainly by electron scattering experiments. 

From the conclusions extracted out of the present calculations, we see 
that the use of deformed Woods-Saxon field in our BCS approach reproduces 
quite well the experimental form factor data in the set of isotopes studied, 
considerably better than the description offered when utilizing a deformed 
harmonic oscillator basis. The simultaneous consideration of the deformation 
and pairing correlations affect largely the Fermi sea depletion and produces 
appreciable smearing of the nuclear surface. 

The investigation of the infiuence of the quadrupole deformation (param- 
eter P) on the surface textures of the above moderately deformed nuclear 
systems, showed that the active orbit occupancies and, subsequently, the de- 
pletion of the Fermi sea for each isotope, are dependent slightly upon the 
parameter f3. The results found for the Fermi sea depletion are more than 
40% (depending on the specific nuclear system) larger than those emerging 
out of the old experimental data, from which one would infer that our method 
overestimates, somehow, the Fermi sea depletion through stronger deforma- 
tion effects in order to compromise for the correlation effects not considered 
by it. However, recent information (experimental and theoretical) tend to 
support larger depletions (of the order of 20-25 %) which are consistent with 
our present estimations. The remaining discrepancy may be clarified from 
specific future experiments of greater precision and larger momentum trans- 
fers. 

The authors would Uke to thank Professor 1. Sick for providing us with the 
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project. 
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Figure 1: Calculated and experimental fH^l^lEl charge form factors for 
the deformed nuclei ^^Mg, HUg, HSi and f^S. The results labeled WS+BCS 
(solid line) include pairing correlations and deformation effects while those 
labeled WS (dashed line) originate from a simple Slater determinant con- 
structed out of a deformed Woods-Saxon basis wave functions. The dashed 
dotted curve, DHO+BCS, comes from a deformed harmonic oscillator field 
employed in the BCS ground state. 
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Figure 2: Charge form factor of ^^Mg (up) and ^^S (down) for various values 
of the deformation parameter /? entering the deformed WS basis. Pairing 
correlations are also taken into consideration through the BCS ground state. 
The values of (3 comes from Ref . [2^] • 
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Figure 3: Proton occupation probabilities near the Fermi surface for the 
nuclei HMg, ^2Mg, f|Si and fgS. The diffuseness of the Fermi surface reflects 
the pairing correlations and deformation effects. The notation is the same 
with that of Fig. 1. 
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Figure 4: The diffused Fermi surface and tlie Fermi sea depletion for i2^S 
(up) and ^^S (down) produced by varying tlie deformation parameter /? in 
a deformed Woods-Saxon basis (see Table 1). The strength for proton-pair 
(neutron-pair) interaction in the BCS procedure is fixed at the value App — 
Ann = 12.84/^^2 MeV. 
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